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Abstract 

It is investigated the statistical properties of random walks evolving on real configurations of 
a crumpled wire rigidly jammed in two dimensions. These crumpled hierarchical structures with 
complex topology are obtained from a metallic wire injected at a constant rate into a transparent 
planar cell of 20cm of diameter. The observed diffusion is anomalous with an exponent very close 
to that obtained at the threshold of two dimensional percolation. A comparison of the system 
studied in this paper with other systems of physical interest is also made, and an experimental 
consequence of our results is discussed. 
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I. INTRODUCTION 

The study of crumpled structures is a topic of increasing scientific interest in physics 
and correlated disciplines Q], as well as in the mathematical and technological domains 
In the last years, theoretical and experimental aspects of the condensed matter physics of 
crumpled sheets have been investigated in many areas of study, e.g. acoustic emission jsj], 
continuous mechanics (|, growth models [3], packing problems ja], polymer, membrane, and 
interface physics Q], universality jiol ]. among others. However, crumpled structures with 
different topologies, as exemplified by a squeezed ball of wire, have been much less studied 
in the physics literature. Geometrical, statistical, and physical aspects of crumpled wires 
( CW) in 3d space were examined fifteen years ago from the point of view of experimental work 
and analogue simulations, and in particular some robust scaling laws and fractal dimensions 
associated with these disordered systems were reported [llj . More recently, it was discovered 
that the irreversible process of injection of a long metallic wire in a two dimensional cavity 
lead to crumpled structures exhibiting several scaling phenomena Q). In the present work, 
we report results of an extensive analysis of diffusion on these 2d crumpled structures of 
wire obtained by irreversible squeezing of macroscopic pieces of copper wires within a two- 
dimensional transparent cell. 

The outline of this article is as follows. In Sec. |H]we describe the essential details from 
the experiment and from the simulations, and in Sec. IHII we present and discuss our results. 
In Sec. II VI we summarize our major conclusions. 



II. EXPERIMENTAL DETAILS AND SIMULATIONS 

The experimental apparatus used in our work to obtain the configurations of 2d CW 
consists of a transparent cell formed by the superposition of two discs of plexiglass with a 
total height of 1.8cm, an external diameter of 30cm and an internal circular cavity with 
R = 10.0cm of radius and 0.11cm of height, which can accommodate configurations of a 
single layer of CW of 0.10cm of diameter. The cavity of the cell was polished and the 
copper wire used in the experiments (#19AWG, with diameter ( = 0.10cm) had a varnished 
surface, in order to reduce the friction. Cavity and wire operated in dry regime, free of any 
lubricant. Two radial channels were made to provide the injection of wire into the cell at 
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FIG. 1: Photograph of the 2d injection cell used in the experiments discussed in this paper. 
The internal cavity contains a typical CW configuration (in this particular case for the maximum 
packing fraction, with 14.0% of the area of the cavity occupied by the CW). All configurations of 
CWin 2d present a hierarchical cascade of loops with several types of contact wire-wire leading to 
a complex topology. See Sec. [0] for detail. 

the angle of 180° , as indicated in Fig^ Digital photographs of the CW configurations were 
taken with an Olympus C-3040ZOOM digital camera with resolution of 3 megapixels, which 
was connected with a PC, and assembled 30cm over the cell. To avoid picture artifacts by 
light reflections a cylindrical paper screen was placed around the cell. 

In order for the reader to develop insight about the nature of the crumpling processes 
that we are considering in this paper, we show in Fig0]a typical 2d configuration of CW at 
the limit of maximum packing corresponding to an occupation of 14.0% of the total area of 
the cell. The packing fraction is given p = where L is the length of the wire introduced 
in the cavity. The maximum length of wire which is possible to introduce in this case is 
L m ax = 440cm (ensemble average), corresponding to p max = 0.14. In fact, this maximum 
packing fraction does not vary appreciably, because it is controlled by the (approximately 

universal) ratio ( shear modulus j Q ^ ^ ne w j re The maximum packing fraction is 

' V modulus of elasticity/ !-Jfl J i — i ° 

also independent of the degree of lubrication of the cavity |l2|. Each injection experiment 
begins fitting a straight wire in the opposite channels and subsequently pushing manually and 
uniformly the wire on both sides of the cell toward the interior of the cavity. The injection 
velocity at each channel in these experiments was of the order of lcm/ s. For occupations 
close to Pmax, the crumpled structures become rigid, the difficulty in the injection increases, 
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and the injection velocity goes to zero abruptly. However, the observed phenomena are 
widely independent of the injection speed for all interval of injection velocity compatible 
with a manual process. The complex crumpled patterns of wire observed within the cavity 
are basically due to the formation of a cascade of loops of decreasing size. During the 
progressive injection of wire into the cell, the cascade of loops evolves in such a way that 
it is common to observe localized or global rearrangements of the loops previously formed, 
provided the packing fraction is not close to p m ax- In principle, there is a concentration of 
stress and energy in the regions of the loops with smaller curvature radii, and a negligible 
amount of stretching is involved. The reader can also observe that the sharp creases and 
ridges found in crumpling of sheets are absent in the 2d CW shown in Fig. ^ For additional 
details on the preparation and geometric properties of 2d CWwe refer the interested reader 
to reference |12| . 

Using high resolution digital images from an ensemble of five independent and equivalent 
configurations of CW at the limit of maximum packing (Fig. ^) we simulate random walks 
by considering the following rules: (i) The random walker has size ( identical to the width of 
the wire, (ii) Each walk begins in a point chosen at random on the CW. (iii) In regions free 
of contact wire-wire, the random walker at position had equal probabilities to go ahead 
(to Ti + i) or to go back (to rj_i) in unit steps, as in a typical one- dimensional trail, (iv) On 
the other hand, in regions with contact wire-wire, the more complex local topology must 
be considered and the walker can move with equal probability from to one of its closest 
neighbors sites {5} (i.e. a number of neighbors > 2) localized at r^g. For illustration, we 
show in Fig. 121 a detail of a region taken from Fig. ^ in which a large number of contacts 
wire-wire occur. In these regions the walker can jump from a loop to another contiguous 
loop, and this possibility modifies drastically the topology of the random walk. From Fig. |21 
one can notice contacts of variable extent involving e.g. 2, 3, 4, and 5 trails. 

III. RESULTS AND DISCUSSION 

Due to the possibility of diffusion transverse to the trails in the regions with a variable 
number of contacts between loops, an extensive study of the statistical properties of the 
random walks on CWis necessary. Using the rules defined in the end of SeclHl we calculated 
the mean square displacement of a random walker, < r 2 (N) >, as a function of the number 
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FIG. 2: A detail of Fig. ^ showing the geometry of the contacts wire-wire. The distribution of 
points within the trail (wire backbone) represents the positions of a random walk of 10,000 steps. 
It can be noticed that the walk is trapped in a very restricted region of the CW - an indication 
that the process is diffusive, with an anomalous diffusion exponent d w > 2. See Sees. HT1 and IIIII for 
detail. 

of steps N as shown in the log — log plot of Fig. 01 This ensemble average was obtained from 
1,150 walks of 10,000 steps for all configurations of CW considered. We see from the main 
plot of Fig. El that < r 2 (N) > is well described by a power law function: < r 2 (N) >~ N a , 
with the slope a assuming the value 0.66±0.01 (Obtained from the interval 100 < N < 3, 000; 
i.e. for relatively not too large number of steps ["short-time" limit] to avoid the finite size 
effect of the cavity associated with the longer walks.). Thus, the random walk exponent 
Sw WS> = - assumes the anomalous value 3.03 ± 0.05, sensibly larger than the brownian 
exponent d w = 2 for normal diffusion. In order to display clearly the power law best fit 
leading to the determination of the exponent a (continuous line) and the typical fluctuation 
bars in the main plot, we exhibit in this figure only 26 data points. The complete overcrowded 
simulation data is represented in the lower right inset. For the sake of completeness we show 
in the upper left inset of Fig. 01 the dependence of the exponent a with the inverse of the 
number of random walks {j^rw) used in the simulations, for some values of j^rw in the 
interval [150, 1150]. The linear fit (continuous line) in this inset suggests that a converges 
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FIG. 3: The main plot shows the log — log dependence of the mean square displacement for 
random walks on a CW, < r 2 (N) function of the number of steps N (only 26 selected 

points are represented to put in evidence the power law fit and the fluctuation bars. The data 
follows the scaling < r 2 (N) >~ N a , with a assuming the value 0.66 db 0.01 (obtained from the 
interval 100 < N < 3,000). The random walk exponent d^ W ^ = ^ assumes the anomalous value 
3.03 ± 0.05. The upper inset shows the dependence of the exponent a with the number of random 
walks used in the ensemble. The lower inset shows the same variables of the main plot, but for the 
complete simulation data. See Sec. Illll for detail. 

to the value 0.655 within uncertainties of 0.010 at the limit j^rw — ► oo. Obviously, a study 



of the finite-size scaling for the random walks along the line suggested in reference |14 | 
would be interesting in this problem. However, the realization of experiments with cavities 
of several different sizes is presently connected with severe technical difficulties. Figure |4(a)| 



and 4(b)| give for completeness additional information on the sample to sample fluctuations 



on < r 2 (N) >: in Fig. |4(a)| differently from Fig. |4(bjj the saturation on < r 2 (N) > is already 
evident for N close to 10, 000 — 13, 000 steps. We notice that the exponents d w obtained 
from the single CW configurations exhibited in these figures are close to the exponent found 
from Fig. El In particular, from Fig. |4(b)| we obtain d w = 2.86 ± 0.04, for 10 < iV < 10 4 . 
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FIG. 4: (a) and (b) show the mean square displacement of a random walker for two individual 
samples (see text, Section HTfl for detail). 



The value obtained from Fig. |3] for the exponent <% is close to the value 
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(PC) 

dw = 2.87 ± 0.02 obtained for anomalous diffusion on the self-similar incipient in- 
finite percolation cluster (PC) at d = 2 jl^. Although 2d CWand 2d percolation clusters 
have different geometric aspects, both structures possess some similarities: CWhas fractal 
dimension D( cw ^ = 1.9±0.1 jl^ (The reader can observe that this is the experimental value 
obtained from the CW configurations studied in i.e. exactly the same CW configura- 
tions used in the present work.), and PC has a close fractal dimension D^ PC ^ = 1.89 ± 0.03; 
they have structures with the same topological dimension, d? = 1; they have a hierarchy of 
holes or empty spaces without any characteristic size, and they equably can support stress. 

The exponent dl^ W ^ reported here for CWis in agreement with the conjecture of Aharony 
and Stauffer [16J which predicts that d w = D + l, for a system with fractal dimension D < 2. 
According to this conjecture, the observed experimental value D^ cw > = 1.9 ±0.1 \12\ . points 
to a value d^w = 2.9 ± 0.1, which is equal to d!^ W ^ = 3.03 ± 0.05 within the statistical 
fluctuations. 

As an additional insight, it is interesting to compare the exponent d^ w ^ with the diffusion 
exponent for a 2d self-avoiding random walk (SAKW) with a (random) distribution of local 
bridges without mass connecting any two nearest-neighbor sites visited by the SARW. This 
structure which could be speculated in principle as a valid zero-order approximation for 2d 
CW allows both diffusion along the backbone of the SARW as well as diffusion transverse 
to the SARW backbone (via jumps of the walker along the bridges) was studied by Manna 
and Roy ^7]]. This SARW with, bridges has the same fractal dimension of the SARW, that 
is D = 1.33. Manna and Roy found that d w = 2.56 ± 0.01 for the SARW with bridges, a 
value somewhat different albeit not too distant from 3.03 ± 0.05 observed for 2d CW. 

An experimental consequence of the diffusion exponent reported in this paper is connected 
with the dependence of the electrical resistance 9ft of a CWwith its size (radius) R. Here we 
are assuming a CW allowing for conduction along the wire backbone as well as transversal to 
it (in the regions of contact wire-wire). In general, the electrical resistance of an Euclidean 
(^-dimensional system of size R scales as 9ft ~ j^ ea — ft(dLi) — P ^ • For a fractal system, 
the corresponding relation is 3ft ~ jQdw-D)^ w here ]j j s the fractal dimension D < d, and 
d w is the random walk exponent obtained from random walks defined on the fractal [18] . 
The low connectivity of fractals lead to d w > 2, i.e. to an exponent larger than the normal 
(brownian) value. Using the observed values D^ cws> and dw for CW, we predict that the 
electrical resistance for such systems scales as 9ft( cw ) ~ R^ 1A±0A \ whereas normal diffusion 
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for a 2d Euclidean system implies that 3ft ~ R( 2 ~ 2 ^ = R°. Thus, the electrical resistance of 
a CW of the type considered above behaves very close to the electrical resistance of a Id 
resistor (3ft ~ size). For comparison, the expected dependence electrical resistance versus size 
for a percolation cluster is 3ft ~ ^>o.98±o.o5^ obtained from the corresponding values of d!^ ^ 
and D^ PC \ On the other hand, the electrical resistance for a CW allowing for conduction 
only along the backbone of wire (i.e. excluding jumps of current in the contacts wire-wire) 
scales as 3ft ~ (length of wire within cavity) ~ R D ~ i? L9±ai . 

Recently (l9j, an investigation in the fluctuation of the density of states in a quantum 
Heisenberg disordered chain did show its equivalence with the noise spectral density within 



the calculation precision. This result strength the conjecture 20j, |2l( which suggest the 



equivalence between 

from the noise j^JI^I, it would be important to use this conjecture in CW systems to obtain 



_ Doth densities. Moreover, since we can predict the diffusion exponents 
2] , it would be impor 
those exponents in an independent way. 



IV. CONCLUSIONS 

Using extensive numerical simulations we have studied the statistical behavior of a ran- 
dom walker evolving directly on the experimental complex configurations of crumpled wires 
confined in a two dimensional cavity, at the limit of maximum packing fraction (Fig. Q). In 
particular, we have obtained the dependence of the averaged square distance < r 2 (iV) > of 
the random walks as a function of the number of steps N. It is found that these quantities 
scale as < r 2 (N) >~ N a , where a — -j- = 0.66 ±0.01, corresponding to a diffusion exponent 
d^ W ^ = 3.03 ± 0.05, when the packing fraction of wire in the cavity approaches the limit 
value 0.14±0.01. This value for the diffusion exponent is close to that obtained for diffusion 
on the two dimensional incipient percolation cluster [15;]. The exponent dl^ w ^ implies that 
the electrical resistance of crumpled wires scales with the radius (size) as 
a result that can be tested with precise electrical measurements. 
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